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Using neutron pair distribution function (PDF) analysis over the temperature range from 1000 K 
to 15 K, we demonstrate the existence of local polarization and the formation of medium-range, 
polar nanoregions (PNRs) with local rhombohedral order in a prototypical relaxor ferroelectric 
Pb(Mgi/3Nb2/3)03. We estimate the volume fraction of the PNRs as a function of temperature 
and show that this fraction steadily increases from % to a maximum of ~ 30% as the temperature 
decreases from 650 K to 15 K. Below T~200 K the PNRs start to overlap as their volume fraction 
reaches the percolation threshold. We propose that percolating PNRs and their concomitant overlap 
play a significant role in the relaxor behavior of Pb(Mg 1 / 3 Nb2/3)03. 

PACS numbers: 77.84.Dy, 61.43.Gt, 64.70.Kb, 61.12.Ld 



X 



Relaxor ferroelectrics, such as lead magnesium niobate 
Pb(Mg 1/3 Nb 2/3 )03 (PMN), have been widely applied 
due to their high, relatively temperature- independent, 
permittivities. Relaxors are characterized by a frequency 
dispersion in their maximum permittivity temperature 
(Tm ~285 K in PMN), with no macroscopic phase tran- 
sition into a ferroelectric state at Tm These be- 
haviors are fundamentally different from those of nor- 
mal ferroelectrics and are similar to spin-glasses 0. It 
is now believed that local polarization, resulting from 
structural and chemical disorder, plays a crucial role in 
relaxor behavior 0, Q. Several models have been pro- 
posed to explain the nature of the interactions between 
these local polarizations and the mechanism for relaxor 
behavior EIJ51, 0, Hi • These models include superpara- 
electric @, 01 an d dipole glasses 0, as well as random 
field interactions Q ■ The microscopic mechanism behind 
relaxor ferroelectricity is, however, not fully understood, 
and is still controversial due to its high degree of com- 
plexity. 

In this Letter, we report the temperature evolution of 
the local and medium-range crystal structure of PMN 
from 1000 K to 15 K using neutron pair distribution 
function (PDF) analysis. We present evidence for both 
local atomic displacements (local polarization) and for 
medium-range (~ 5 A - 50 A) ordering, called polar 
nanoregions (PNRs). These medium-range correlations 
are modeled using rhombohedral symmetry, enabling for 
the first time an estimate of the temperature dependence 
of the volume fraction of the PNRs. This fraction in- 



creases with decreasing temperature, reaching the three 
dimensional percolation threshold below T^200 K, where 
the PNRs start to overlap. We discuss the implications 
of these findings on the dielectric properties of PMN. 

Evidence for the existence of local polarization in 
PMN below Td ~ 620 K, known as the Burns temper- 
ature jjL was deduced from optical and strain measure- 
ments 0, 0, . The atomic nature of this polarization 
in PMN was also studied using extended x-ray absorption 
fine structure (EXAFS) [H 

[l2j and atomic pair distri- 
bution function (PDF) analysis HQ. Nb K-edge EX- 
AFS measurements on PMN suggested a displacement of 
Nb ions with respect to their oxygen octahedra by about 
0.1 A from room temperature down to 4.5 K. Pb K-edge 
EXAFS measurements indicated strong static disorder 
on the Pb sites, but provided no detailed information 
about static Pb displacements ^lj. This is in contrast 
to neutron PDF measurements on PMN which clearly 
showed that the Pb ions are off-center with respect to the 
O12 cage by as much as 0.5 A [l^l ll^j. In these earlier 
PDF measurements, however, a complete picture of local 
structure in PMN was not obtained due to the limited 
momentum transfer Q and temperature range available. 
Thus, the nature of PNRs and their volume fraction were 
not determined. 

Our measurements were performed on the NPDF 
instrument at the Los Alamos Neutron Science Cen- 
ter (LANSCE). Powder diffraction patterns were cor- 
rected for background, absorption and multiple scatter- 
ing, and normalized using a vanadium spectrum to ob- 
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tain the total scatterin g st ructure function S(Q), using 
the PDFgetN program jl^. The PDF G(r) is obtained 
from S(Q) via the Fourier transform shown in Eq. 1, 

2 rQ m „ 

Gir) = 4irr{p(r) - p ) = - / Q[S(Q) - l]sinQrdQ, 

n Jo 

(1) 

where p(r) and po are the atomic number and average 
number densities, respectively. Since S(Q) includes both 
Bragg and diffuse scattering, the resultant PDF pro- 
vides short-, medium-, and long-range structural infor- 
mation . This technique has been used to study local 
atomic structures |l7| and correlated atomic motions of 
atom pairs |l8| in many materials. 
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FIG. 1: Temperature evolution of PDF spectra of PMN 
within a unit cell, (a) High temperature : T=1000 K, 850 K, 
750 K, 700 K, 650 K, 600 K, 550 K, 500 K, 400 K, and 
313 K, with Q m ax = 25 A -1 . The inset shows a comparison 
between an experimental PDF (symbol) and a model PDF 
(line) at 650 K. See text for details, (b) Low temperature 
: T=250 K, 220 K, 190 K, 150 K, 100 K, 50 K, and 15 K, 
with Q m ax = 30 A -1 . The tick marks indicate (from left to 
right, respectively) Mg/Nb-O, (Pb-O, O-O), Mg/Nb-Pb, and 
(Pb-Pb, O-O, Mg/Nb-Mg/Nb) bond lengths in an ideal cu- 
bic perovskite structure of PMN. The inset shows the cubic 
crystal structure with the Pb ion at body center. 

We first examine the temperature evolution of the low- 
r region of the PDF spectra. Fig. 1 (a) shows spectra 
from 1000 K to 313 K. The tick marks indicate the PDF 
peak positions for an ideal cubic perovskite structure. 
The first tick mark represents Mg/Nb-O, the second Pb- 
O, O-O, the third Mg/Nb-Pb, and the fourth Pb-Pb, O- 
O, Mg/Nb-Mg/Nb bond lengths. At high temperatures 
note that the peaks are quite broad and the Pb-Mg/Nb 
bond is not well defined due to an overlap with the Pb-0 



bond. As the temperature decreases to room tempera- 
ture, the Pb-0 peak gradually resolves into two peaks 
at ~ 2.86 A and 2.5 A. This splitting indicates a dis- 
placement of Pb ions in their O12 cages, which breaks 
12- fold degeneracy of Pb-0 bond distance (see the inset 
of Fig 1(b)). 

In order to understand the origin of these phenom- 
ena we compared the experimental PDF at 650 K with a 
model PDF at the same temperature (shown in the inset 
of Fig. 1 (a)). In the model PDF we assumed no static 
displacements and calculated the widths using thermal 
parameters for Pb, Mg/Nb and O ions obtained from 
our Rietveld refinement at 650 K. Since the area of the 
Pb-0 peak should be conserved, the comparison shown 
in the inset to Fig. 1 indicates that the Pb-0 peak is 
much broader in the experiment than in the model, im- 
plying that significant static Pb local displacement ex- 
ists already at 650 K. Thus, the gradual splitting of the 
Pb-0 peak at lower temperatures is mainly due to a de- 
crease of thermal vibrations. Static Pb displacements in 
PMN above the Burns temperature have also been de- 
duced from Raman measurements [l9| . 

Below room temperature (Fig. 1 (b)) the experimen- 
tal PDFs reveal more interesting structure. First, the 
Mg/Nb-O peak around 2.02 A splits into two peaks, sug- 
gesting local Mg/Nb displacements along the (111) direc- 
tion. Note that this peak shows almost no temperature 
dependence between 250 K and 15 K, consistent with 
Nb K-edge EXAFS measurements 0. In general, PDF 
peaks become sharper with decreasing temperature due 
to decreasing thermal motion. Thus, this unusual tem- 
perature dependence indicates an increasing distribution 
of static Mg/Nb displacements with decreasing temper- 
ature. Second, the Pb-0 peak around 2.85 A splits into 
three peaks at ~ 2.45 A, 2.85 A, and 3.33 A. Third, the 
Pb-Mg/Nb peak around 3.5 A splits into two Pb-Mg/Nb 
bonds at ~ 3.33 A and 3.64 A mostly due to the dis- 
placement of Pb ions. Thus, the peak around 3.3 A has 
both Pb-0 and Pb-Mg/Nb components. The splitting of 
the Pb-0 and Pb-Mg/Nb peaks below room temperature 
suggests that the Pb ions may locally shift mostly along 
(100) or (111) directions. Interestingly, the peak around 
4.05 A (lattice constant of PMN) is well defined; thus, 
the overall cubic symmetry is maintained although the 
local cubic symmetry is not. 

We now address the medium-range correlations 
(PNRs) between these local atomic displacements. 
Fig. 2 (a) shows the PDF spectra from 1000 K to 15 K in 
the r-range up to 16 A, and Fig. 2 (b) shows the temper- 
ature dependence of the peak intensities in the doublet 
between 8.0 A < r <9.8 A. Here both peaks have con- 
tributions from almost all possible ionic pairs. The solid 
line is the calculated temperature dependence of the peak 
height using the average Pb, Mg/Nb, and O thermal pa- 
rameters from Bonncau et al. |20j , which are reproduced 
in the inset. Note the strong deviation of the data below 
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FIG. 2: (a) Temperature dependence of the PDF spectra of 
PMN at temperatures T=1000 K, 750 K, 650 K, 500 K, 313 K, 
250 K, 190 K, 100 K, 15 K. All PDFs were obtained with 
Qmax=25A -1 . (b) PDF peak height in the doublet between 
8.0 A < r <9.8 A as a function of temperature. The line 
is the expected average PDF peak height using the thermal 
parameters in the inset. The inset shows the thermal parame- 
ters of Pb (square), Mg/Nb (circle) and O (triangle) [2(1 . (c) 
Comparison of rhombohedral (R-model) PDF with the exper- 
imental PDF at 650 K. R and Exp represent the R-model and 
experimental (paraelectric cubic) PDFs, respectively. 



T~250 K, an indication of PNRs, as we now show. 

We modeled the PNRs using rhombohedral symmetry 
(R-model: space group R3m), wherein the Pb, Mg/Nb 
and oxygen octahedra are displaced along the (111) direc- 
tions, assuming rigid oxygen octahedra. Similar models 
of rhombohedral correlations with (111) atomic displace- 
ments 0, H3| and non-collinear (100) displacements 
that average to rhombohedral symmetry [24| have been 
proposed in Pb(Mg 1 / 3 Nb 2 /3)03 and Pb(Sc 1 / 2 Ta 1 /2)03, 
respectively. Attempts to model our data using displace- 
ments of orthorhombic symmetry (with the atoms and 
oxygen octahedra displaced along (110) directions with 
the same magnitudes as in the rhombohedral model) were 
unsuccessful, with some PDF peaks being too strong and 
others too weak. 

In Fig. 2 (c) we compare the R-model PDF with the 
experimental PDF at 650 K. For the R-model PDF calcu- 
lation the thermal parameters were taken from Bonncau 
et al. |20| . and the following atomic positions were used, 
determined to best simulate the differences between low 
temperature PDFs and the PDF at 650 K: Pb (-0.0392, 
-0.0392, -0.0392), Mg/Nb (0.5062, 0.5062, 0.5062) and O 
(0.5308, 0.5308, 0.0308). 



The PDF peaks of the R-phase are quite distinct from 
those of the paraelectric cubic phase. In fact, the relative 
intensities of the PDF peaks in the R-phase and in the 
paraelectric cubic phase resemble the low-temperature 
and high-temperature end members shown in Fig. 2 (a), 
respectively. Therefore, if we assume that PNRs are dis- 
persed as "islands" in the paraelectric cubic lattice |2lj . 
growing in size and volume fraction with decreasing tem- 
perature, we expect that features of the R-model PDF 
will become more evident with decreasing temperature. 
This is exactly what wc observe in the experimental 
PDFs shown in Fig. 2 (a). 



(b) T=250K, R3m:20% 
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FIG. 3: Comparison of the "two-phase" model PDF 
with the experimental PDFs at (a) T=500 K, (b) 
T=250 K, (c) T=15 K. The model PDF is calculated from 
G m = a(T)G R +(l-Q(T))G c , where G R and G c are model 
PDFs for the rhombohedral and paraelectric phases. a(T) is 
the volume fraction of the rhombohedral phase, (d) The vol- 
ume fraction increases with decreasing temperature, reaching 
30% at 15K. The inset shows the correlation length of local 
polarization as a function of the temperature (taken from Xu 
et al. [H). 

The volume fraction of the PNRs was estimated us- 
ing a simple model PDF: G m = a(T)G R +(l-a(T))G c , 
where Gr and Gc are rhombohedral and paraelectric 
cubic model PDFs, respectively, and a(T) is the vol- 
ume fraction of the rhombohedral phase as a function 
of temperature. For the R-model PDF calculations the 
atomic coordinates given above were used for all temper- 
atures. For the paraelectric cubic phase, we first intro- 
duced static displacements of ions in a 10 x 10 x 10 unit 
cell via a reverse Monte-Carlo (RMC) fit to the the exper- 
imental PDF at 650 K, using the DISCUS program [26j . 
These static displacements were then fixed, and the cu- 
bic model PDF was calculated at various temperatures 
using the thermal parameters from Bonneau et al. |2Cj . 
(The same thermal parameters are used for both cubic 
and rhombohedral phases.) After calculating both para- 
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electric cubic and rhombohedral PDFs at a given tem- 
perature, a(T) was adjusted to obtain the best match 
to the corresponding experimental PDF. Figures 3 (a-c) 
show these comparisons at 500 K, 250 K, and 50 K in the 
r-range up to r=16 A. Considering the simplicity of the 
model PDF calculations, the overall agreement at various 
temperatures is very good. 

In Fig. 3(d) we show the temperature dependence of 
a(T), along with the correlation length of local polar- 
ization determined by Xu et al. using neutron elastic 
diffuse scattering |2f|, shown in the inset. The volume 
fraction steadily increases with decreasing temperature, 
reaching ~ 30% at 15 K. The correlation length £ of local 
polarization increases from £ ~ 15 A to £ ~ 60 A. Note 
that around T~200 K, where the correlation length dras- 
tically increases, a(T) reaches the percolation threshold 
P c ~28% for spherical objects in three dimensions [27j |. 
This implies that the PNRs start to overlap with each 
other below T~200 K. This correlates well with the 
abrupt increase in the correlation length of local polar- 
ization to £ ~ 50 - 60 A, determined independently. 

Neutron scattering measurements by Wakimoto et al. 
[2^ suggest that locally, each PNR maintains a stable 
spontaneous polarization below T~220 K where the cor- 
relation length is about ~ 50-60 A. Further growth of 
these large PNRs (beyond ~ 60 A) may be hindered 
by electric and elastic energy barriers, which are pro- 
portional to the size of the PNR. In an applied electric 
field, however, the polarizations in the separate PNRs 
are aligned, and the separate PNRs form a macro polar 
domain. Once the macro polar domain is formed, the 
dipolcs in the domain do not randomly orient even after 
the electric field is turned off. 

This picture helps to explain some basic behavior 
found in PMN, such as hysteresis in the field-induced 
polarization [2jj, the field- induced rhombohedral phase 
transition below 220 K j^l and the anomaly in the di- 
electric permittivity near 212 K upon field heating after 
zero-field cooling [Jj. Our measurements, therefore, pro- 
vide direct local structural evidence for the importance 
of percolating PNRs in explaining many features of the 
relaxor behavior found in PMN below the Burns temper- 
ature. 
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